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I. 


INTRODUCTION 


A.  Background 

In  1959  K.  Stewartson1  published  a  theory  on  the  motion  of  a  top 
containing  a  rigidly-spinning  inviscid  liquid  in  a  cylindrical  cavity; 
this  theory  is  equally  applicable  to  shell  and  gyroscopes.  In 
conjunction  with  this  theory  the  gyroscope  provides  an  excellent 
laboratory  tool  for  studying  the  yawing  motion  of  a  spinning  liquid- 
filled  shell. 

Stewartson  predicted  the  existence  of  instabilities  in  the  motion 
of  the  top  due  to  resonances  between  oscillations  of  the  liquid  and 
nutation  of  the  top.  Ward2  experimentally  confirmed  the  presence  of 
two  of  the  predicted  instabilities  in  a  gyroscope,  though  his  results 
differed  in  detail  from  the  prediction.  Experiments  conducted  at  the 
BRL  by  Karpov3**4  with  liquid-filled  shell  and  with  liquid-filled  gyro¬ 
scopes  also  confirmed  the  presence  of  instabilities,  but  showed  a 
strong  dependence  of  the  growth  rate  of  the  instabilities  on  the 
viscosity  of  the  liquid.  Wedemeyer5 6*3  developed  a  theory  of  viscous 
corrections  to  Stewartson1 s  theory  which  predicted  viscous  effects  of 
the  nature  of  those  found  by  Karpov;  his  comparison  of  the  theory  with 
Karpov's  experiments  was  flawed,  as  discussed  below.  Frasier  and 


1.  K.  Stewartson "On  the  Stability  of  a  Spinning  Top  Containing 
Liquid .  "  Journal  of  Fluid  Mechanics ,  Vol.  5.  Part  4 ,  September 
1959 3  pp.  577-592. 

2.  G.  N.  Wardj  Appendix  to  Reference  1. 

3.  B.  G.  Karpov 3  "Dynamics  of  a  Liquid-Filled  Shell:  Resonance  and  the 
Effects  of  Viscosity s"  BRL  Report  No.  12793  U.S.  Army  Ballistic 
Research  Laboratory 3  Aberdeen  Proving  Ground 3  Maryland 3  May  1965. 

AD  468654. 

4.  B.  G.  Karpov 3  "Liquid  Filled  Gyroscope:  The  Effect  of  Reynolds 
Number  on  Resonance 3"  BRL  Report  No.  1302s  U.S.  Army  Ballistic 
Research  Laboratory ,  Aberdeen  Proving  Ground.  Maryland.  October 
1965.  AD  479430. 

5.  E.  H.  Wedemeyer 3  "Dynamics  of  Liquid-Filled  Shell:  Theory  of 
Viscous  Corrections  to  Stewartson1 s  Stability  Problem 3"  BRL  Report 
No.  1287 j  U.S.  Army  Ballistic  Research  Laboratory Aberdeen  Proving 
Ground^  Maryland3  June  1965.  AD  472474. 

6.  E.  H.  Wedemeyer  3  " Viscous  Corrections  to  Stewartson* s  Stability 
Criterions"  BRL  Report  No.  13253  U.S.  Army  Ballistic  Research 
Laboratory s  Aberdeen  Proving  Grounds  Maryland 3  June  1966. 

AD  489687. 
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Scott7  developed  a  theory  analogous  to  Stewartsonfs  for  a  gyroscope 
with  a  solid  rod  concentric  with  the  cylindrical  cavity.  They  predicted 
that  the  presence  of  the  rod  primarily  would  change  the  natural  frequen¬ 
cies  of  the  liquid,  the  motion  itself  remaining  similar  to  that  without 
the  rod.  Frasier8  applied  Wedemeyer fs  theory  of  viscous  corrections  to 
this  case  and  performed  gyroscope  experiments  to  check  the  theory. 

Karpov  and  Wedemeyer  used  the  works  mentioned  above,  plus  others, 
as  the  foundation  for  Ref.  9,  the  handbook  for  liquid-filled  projectile 
design.  Ref.  9  presents  theoretical  and  experimental  results  as  well 
as  operational  formulas  intended  for  use  by  designers  in  assessing  the 
stability  of  liquid-filled  projectiles.  Since  the  preparation  of  the 
handbook,  D'Amico10,  Scott  and  D'Amico11,  and  Kitchens  (unpublished) 
have  performed  experiments  in  which  growth  rates  of  the  unstable  motions 
of  liquid-filled  gyroscopes  were  measured.  The  object  of  most  of  these 
experiments  was  to  observe  changes  in  growth  rate  with  large  amplitude 
of  the  motion,  and  no  comparison  has  been  reported  between  the  results 
of  these  experiments  and  the  linear,  small  amplitude  theory  of 
Stewartson-Wedemeyer ,  though  such  comparison  is  possible. 

B .  Purpose  of  This  Report 

Wedemeyer 's  theory  of  viscous  corrections  for  solid-body  rotation 
consists  of  two  distinct  parts:  (1)  a  theory  of  viscous  effects  on 
oscillations  of  the  liquid,  and  (2)  a  theory  of  the  response  of  the 
gyroscope  to  the  oscillations.  At  the  time  Karpov4  was  reporting  his 
final  experiments,  only  the  second  part  of  the  theory  had  been 
developed--Wedemeyer  could  predict  the  motion  of  the  gyroscope  given 
the  oscillations  of  the  liquid,  but  could  not  accurately  predict  the 
frequencies  and  damping  factors  of  these  oscillations. 

7.  J.  T.  Frasier  and  W.  E .  Scotty  " Dynamics  of  a  Liquid-Filled  Shell: 
Cylindrical  Cavity  With  a  Central  Rod,11  BRL  Report  No.  1391,  U.S . 
Army  Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground, 
Maryland,  October  1968.  AD  667365. 

8.  J.  T.  Frasier,  irDynamics  of  a  Liquid-Filled  Shell:  Viscous  Effects 
in  a  Cylindrical  Cavity  With  a  Central  Rod,11  BRL  Memorandum  Report 
No.  1959,  U.S.  Army  Ballistic  Research  Laboratory,  Aberdeen 
Proving  Ground,  Maryland,  January  1969.  AD  684344. 

9 .  Engineering  Design  Handbook,  Liquid-Filled  Projectile  Design,  AMC 
Pamphlet  No.  706-165,  U.S.  Army  Materiel  Development  and  Readiness 
Command,  Washington,  D.  C.,  April  1969.  AD  853719. 

10.  W.  P.  D*Amico,  nInertial  Mode  Corrections  for  the  Large  Amplitude 
Motion  of  a  Liquid-Filled  Gyroscope,  11  PhD  Thesis,  University  of 
Delaware,  Newark,  Delaware,  June  1977. 

11.  W.  E.  Scott  and  W.  P.  DfAmico,  11  Amplitude- Dependent  Behavior  of  a 
Liquid  Filled  Gyroscope,  ,r  J.  Fluid  Mech. ,  Vol.  60,  Part  4,  October 
1973,  pp.  751-758. 
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Karpov  used  his  measurements  of  the  motion  of  the  gyroscope  and 
the  second  part  of  the  theory  to  find  the  liquid  frequencies  and  damping 
factors  at  resonance,  as  explained  in  Section  5  of  Ref.  4.  Fig.  3  of 
this  reference  shows  the  resulting  damping  factors.  Karpov  used  these 
values  and  the  second  part  of  the  theory  (Eq.  (5)  in  Ref.  4)  to  compute 
the  growth  rate  of  the  motion  of  the  gyroscope  away  from  resonance. 

This  calculated  growth  rate  is  shown  by  the  solid  curves  in  Figs.  4  and 
5  of  Ref.  4;  the  circles  in  those  figures  are  data.  Because  it  was 
necessary  to  use  the  measured  resonance  growth  rate  itself  in  predicting 
other  growth  rates,  the  theoretical  results  could  only  predict  the  shape 
of  the  growth  rate  (vs.  an  appropriate  parameter)  curve  in  the  neighbor¬ 
hood  of  resonance.  There  was  no  check  of  the  magnitude  of  the  growth 
rate  at  any  frequency. 

Fig.  1  of  Ref.  6  presents  the  same  information  as  Figs.  4  and  5  of 
Ref.  4.  Thus,  despite  the  statements  at  the  bottom  of  p.  22  of  Ref.  6, 
these  figures  present  only  a  very  limited  check  of  theory  with  experi¬ 
ment.  These  same  diagrams  are  presented  as  Fig.  6-2  of  Ref.  9;  in  the 
text  preceeding  the  figure  (Section  6-4)  a  succinct  explanation  is  given 
of  the  way  in  which  the  solid  curves  were  actually  computed.  Note  that 
in  this  explanation,  the  relation* 

[  |  D |  /  (aL)  ] ^  =  1.6  x  10"3 

is  given  as  descriptive  of  the  cylinder.  The  meaning  of  this  quantity 
will  be  explained  later;  what  is  of  significance  here  is  that  for  the 
reported  experiment,  it  was  not  constant,  varying  by  more  than  a  factor 
of  two  over  the  experimental  range.  Thus,  the  first  part  of  Wedemeyer's 
theory  has  not  been  checked  at  all  against  Karpov's  data,  and  the  second 
part  has  not  been  checked  consistently. 

Frasier  modified  the  theory  for  the  rodded  case^  by  Wedemeyer's 
viscous  correction  technique.  The  comparison  between  the  modified 
theory  and  experimental  data  is  flawed,  however,  because  the  theory  was 
not  correctly  evaluated.  A  discussion  of  the  nature  of  the  error  is 
deferred  until  Section  III. 

The  experiments  conducted  after  the  release  of  the  liquid-filled 
projectile  handbook.  Ref.  9,  have  not  been  used  to  check  Wedemeyer's 
theory,  in  part  because  they  were  intended  for  other  purposes.  Thus, 
because  of  the  difficulties  with  Karpov's  and  Frasier's  comparisons, 
the  situation  exists  that  the  theory,  which  forms  a  substantial  part 
of  the  handbook,  has  never  been  adequately  checked  against  experimental 
data,  though  data  are  available. 

The  primary  purpose  of  this  report  is  to  present  a  comparison  of 
the  theory  with  all  of  the  data  from  gyroscope  experiments  conducted 

*  Definitions  are  given  in  List  of  Symbols  section j  p.  33. 
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at  BRL,  so  that  a  proper  assessment  of  its  validity  can  be  made.  In 
addition,  some  significant  misstatements  and  errors  in  equations  in 
Ref.  9  are  pointed  out  and  corrected. 


I I .  THEORY 


A.  Equations  of  Motion  of  the  Gyroscope 

A  brief  sketch  of  the  theory  is  given  here;  the  reader  is  referred 
to  Refs.  1,  6,  7,  and  9  for  details.  The  purpose  here  is  to  show  where 
the  theory  is  susceptible  to  experimental  confirmation. 


The  gyroscope,  shown  schematically  in  Fig.  1,  consists  of  (1)  a 
rotor  spinning  around  its  symmetry  axis  at  angular  speed  ft,  and 
(2)  mounts  which  allow  the  rotor  (z'  )  axis  to  yaw  about  a  fixed  line 
(z-axis)  through  the  "point  of  support”  (the  point  defined  by  the  inter¬ 
section  of  lines  drawn  through  the  gimbal  pivots).  Let  the  (x,y,z) 
coordinate  system  be  an  inertial  system  with  the  x-axis  passing  through 
the  outer  gimbal  pivots  and  the  z-axis  parallel  to  gravity.  Let  the 
non-rotating  ex'  ,/  ,  z'  )  coordinate  system  be  attached  to  the  inner 
gimbal  mount  in  such  a  way  that  when  the  z  and  z'  axes  coincide,  the  x 
and  x'  and  y  and  /  axes  also  coincide. 


Fig.  2  shows  coordinates  (the  x'  and  /  axes  being  omitted  for  the 

sake  of  clarity).  The  z'  -axis  is  coning  about  the  z-axis  at  a  small 

yaw  angle,  aft),  where  a  <<  1 .  The  yawing  motion  of  the  gyroscope  is 

described  by  the  component  angles  0  and  0  of  the  projection  of  the 

x  y 

end  point  of  a  unit  vector  in  the  z'  -direction  upon  the  x-y  plane,  as 

shown  in  the  figure.  We  combine  these  components  in  the  complex  yaw 

0=0  +  i  0  . 

x  y 


The  liquid  is  confined  in  a  cylindrical  cavity  whose  symmetry  axis 
coincides  with  that  of  the  rotor;  the  radius  of  the  cavity  is  a,  and 
the  length  is  2c.  For  a  partially-filled  cavity,  the  fraction  of  fill 

is  1  -  b2/a2,  where  b  is  the  inner  radius  of  a  concentric  annulus  with 
volume  equal  to  that  of  the  liquid.  In  the  case  of  a  rodded  cylinder, 
the  air  core  is  replaced  by  a  rigid  rod  of  radius  d. 

There  are  two  torques  acting  on  the  gyroscope,  the  first  caused  by 
gravity  acting  on  the  mass  of  the  gyroscope,  and  the  second  caused  by 
the  liquid.  Assuming  small  0,  we  denote  the  first  moment  by  Mq0  and 

the  second  moment  by  T0,  where  M  is  a  real  constant  (  =  0  when  the 

°  * 

pivot  point  and  center  of  mass  coincide)  and  F  (0,  d0/dt)  is  a  complex 
function.  The  equation  of  motion  of  the  gyroscope  when  | 0 |  is  small 
(see,  e.g.,  pp.  2-14  and  3-11  in  Ref.  9)  is 
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(1) 


T  d20/dt2  -  iLfi  dG/dt  -  MB  =  -  f  (6,  d0/dt)  0  , 

where  L  is  the  axial  moment  of  inertia  of  the  gyroscope,  and  T  is 
the  transverse  moment  of  inertia  about  the  pivot  point.  The  moments 
of  inertia  about  the  x  -  and  y  -axes  must  be  approximately  equal  here; 
see  p.  2-13  of  Ref.  9  for  discussion. 

* 

Because  of  the  dependence  of  Y  on  0,  Eq.  (1)  is  non-linear;  the 
solution  to  an  initial  value  problem  cannot  in  general  be  constructed 
by  linear  combination  of  independent  solutions.* 

B .  Exponential  Solutions 

For  certain  values  of  the  complex  constant  x,  possible  solutions 
of  Eq .  (1)  are 


0  =  exp  (i  x  Q  t)  .  (2) 

These  solutions  constitute  the  bases  of  the  Stewartson-Wedemeyer 
analyses,  particularly  in  regard  to  establishing  criteria  for  stable 
motion.  The  permissible  values  of  x  are  found  by  substituting  Eq.  (2) 
into  Eq.  (1),  which  gives 

T  t2  -  L  x  +  (M  /ft2)  =  r  (t) /ft2  .  (3) 

Once  the  liquid  moment  function,  V  (x),  is  provided,  Eq.  (3)  is  solved 
for  x. 

Let  x  =  xD  +  i  xT.**  Eq.  (2)  then  represents  an  angular  motion, 

K  1 

here  called  coning  motion,  of  constant  frequency,  xnQ,  which  decays 

with  time  at  the  rate  x^-ft  if  Xj.  is  positive,  but  grows  in  amplitude 

and  becomes  unstable  if  x^  is  negative.  The  quantity  -  x^  is  called 

the  "yaw  growth  raten,  and  its  measurement  is  the  primary  objective 
of  gyroscope  experiments. 


*  In  this  regard  the  statement  following  Eq.  (4-17)  in  Ref.  9  is 
misleading . 

**  In  the  original  inviscid  perturbation  theory  of  Stewartson x  is 
real  for  certain  regions  in  parameter  space ;  see  the  un-numbered 
equation  preceeding  Eq.  (5.12)  in  Ref.  1. 
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The  forcing  term  on  the  right-hand  side  of  Eq.  (1)  is  (1)  zero  if 
the  gyroscope  is  empty,  and  (2)  negligible  if  is  not  near  any 

natural  frequency,  t  ,  of  free  oscillation  of  the  liquid  (see  bottom 

p.  10  of  Ref.  6),  i.e.,  away  from  resonance.  For  either  of  these 
conditions,  Eq.  (3)  has  two  real  solutions,  denoted  by  and  xp  (nuta- 

tional  and  precessional  frequencies,  respectively,  of  the  gyroscope): 

Tn  =  [L/  (2T)  ]  [1  +  (1  -  B)?2]  (4) 

Tp  =  [L/ (2T) ]  [1  -  (1  -  3)^] 

where 

3  =  4  Mq  T/(LQ)2  .  (5) 

In  these  situations,  0 (t)  is  a  linear  combination  of  exp  (i  fit)  and 
exp  (i  Tp  ft  t)  . 

The  basis  for  the  prediction  of  yaw  growth  rate  is  the  premise 
that  near  resonance  the  motion  of  the  gyroscope  with  liquid  is 
dominated,  over  a  significant  time  interval,  by  a  single  mode  of  the 
form 


0=0  exp  (i  td  ft  t)  exp  (-  tt  ft  t)  ,  (6) 

OK  1 

where  ij  is  negative.  This  premise  is  supported  by  experimental 

evidence  and  the  conclusions  of  Stewartson-Wedemeyer  theories.  Then 
the  amplitude  of  the  yaw  should  behave  as 

In  (modulus  [0  ]!)  -  const  -  tt  0.  t  ,  (6a) 

meas  1 


which  represents  a  straight  line  whose  slope  provides  t^. 

C.  Evaluation  of  the  Liquid  Moment 

Determination  of  the  moment  exerted  by  the  liquid  on  the  gyroscope 
requires  solving  the  appropriate  equations  of  motion  for  the  liquid, 
subject  to  boundary  conditions  determined  by  the  position  and  motion 
of  the  rotor.  The  solution  to  these  equations  yield  pressure; 
integration  of  the  pressure,  with  the  appropriate  moment  arm, 
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over  the  inner  surface  of  the  rotor  would  give  the  liquid  moment.*  To 

the  authors'  knowledge  no  one  has  actually  carried  out  this  integration. 

The  integration  can  be  avoided  in  the  region  of  resonance  (t„  x  )  by 

R  oo  J 

employing  an  approximation  of  Stewartson1 . 

C. 1 .  Inviscid  Theory:  Approximation  Near  Resonance. 

The  moment  exerted  by  an  inviscid  liquid  is  very  small  unless 

the  coning  frequency  x  is  near  one  of  the  liquid  natural  frequencies 

x  .  If  x  is  near  x  ,  i.e.,  xn-x  I  and  I  xT  I  «1,  then  T  can  be 
o  o  1  R  o 1  1  1 1 

approximated  in  Eq.  (3)  by 

r  (x)  -  n2  D  (xo)  /  (x  -  xo)  ,  (7) 


where  D  is  the  residue  of  a  Laurent  series  expansion  about  t 
is  convenient  to  define  a  non-dimensional  quantity  R: 


t  .  It 
o 


r2  =  -  D  (J0)  c  /  (p  a6) 


(8) 


The  eigenfrequencies  and  the  residues  R  depend  only  on  c/a,  b/a, 

and  the  azimuthal,  axial,  and  radial  wave  numbers  of  the  liquid  wave: 
m,  j,  and  n.  Only  singly-periodic  waves  (m  =  1)  contribute  to  the 
moment,  and  only  they  are  considered  hereafter. 

Tables  of  tq  and  2R  vs  c/[(2j  +  1)  a]  and  b2/a2  are  presented  in 

Ref.  1**  for  cylinders  without  a  central  rod,  and  in  Ref.  7  (d  replacing 
b***)  for  cylinders  with  a  central  rod.  Expanded  versions  of  both 
tables  are  presented  in  Ref.  9. 

C. 2.  Viscous  Correction  for  Large  Reynolds  Number 

For  quantitative  purposes  the  effect  of  viscosity  cannot  be 
neglected  in  the  prediction  of  eigenfrequencies,  and  consequently,  of 
yaw  growth  rates.  Wedemeyer6  demonstrated  that  for  large  Reynolds 
numbers,  Re  (=  aft2/v),  a  viscous  correction  to  the  inviscid  eigen- 
frequency  can  be  derived.  In  the  formalism  of  his  treatment,  the 


*  We  note  that  minus  signs  should  appear  before  the  integrals  on  the 
right-hand  sides  of  Eqs.  (3-31)  and  (3-30a)  in  Ref .  9. 

**  In  Ref .  13  Tables  2-5_,  the  column  labeled  "i?"  should  be  labeled 
" 2R 


***  The  rrb2/a2tr  labeling  on  pp.  53-67  of  Ref .  7  should  be  replaced  by 
”d2/a2”. 
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viscous  eigenfrequency  (=  tqo  +  i6)  is  complex,  with  a  viscous 

damping  factor  6;  tqo  differs  from  t  .  The  presence  of  boundary  layers 

on  the  cylinder  walls  effectively  changes  the  cavity  dimensions  to  a-6a, 
c-6-t,  and  d+5d,  where  the  complex  increments  6a,  6c,  and  6d  are  given 
by* : 


6a/a  =  (2Re)~li  (1  +  i)  (1  - 


(9a) 


6c/c  = 


(2Re)  2  a 


2  ( 1  -t  )  c 
o 


(1-i)  (3-xo) 

(1+T  )  2 
o' 


(1+i) (1+tq) 
(3-t  )  2 


(9b) 


6d/d  =  a/d  (6a/a) 


(9c) 


The  viscous  eigenfrequency  is  then  obtained  formally  as  a  first- 
order  correction  to  x  .  For  a  cylinder  with  no  central  rod. 


3  T 

C  0 

6a  6c 

l 

2b2 

9  T 

0 

a(2j+l)  3 [c/a(2j+l) ] 

a  c 

'T 

1 

a2 

3(b2/a2)_ 

(10) 


For  a  cylinder  with  a  central  rod 


8  T 


T  =  T  + 
O  V  o 


a(2i+1)  9 (ca_1 

{ 2j+l}-1 ) 

2d2  3  To 

6a  6d 

a2  9(d2/a2) 

a  +  d 

6a  6c 
a  c 


(11) 


*  Wedemeyer  denotes  the  argument  in  the  right-hand  side  of  Eqs .  (17) 
and  (19)  of  Ref .  6  by  x,  calling  it  the  n dimensionless  frequency  of 
oscillation Here  he  is  dealing  with  free  oscillations  of  the 
liquid j  so  t  should  be  set  equal  to  x  .  In  Ref .  93  Eqs .  (6-13)  and 

(6-14)  also  contain  the  argument  x,  but  the  succeeding  statement 
,fwith  x  =  a)/Qj  the  dimensionless  yawing  frequency "  is  incorrect; 
however j  in  the  example  at  the  top  of  p.  6-8>  x  is  properly  set 
equal  to  x  • 
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The  effect  of  viscosity  on  the  liquid  moment  is  then  approximated 
in  Wedemeyer's  theory  by  replacing  tq  in  Eq.  (7)  by  xqv,  leading  to  the 

following  new  approximation  to  Eq.  (3) : 

T  t2  -  L  x  +  (Mo/fi2)  =  D  (to)  /  [x  -  (too  +  i  6)]  .  (12) 


D.  Solutions  for  x 

Eq.  (12)  has  three  complex  roots.  Because  of  the  approximation, 
Eq.  (7),  made  in  evaluating  the  liquid  moment,  only  those  roots  with 
XR  ^  Too  are  aPPlicable-  If  t  **  or  x  »  x-n,  the  left  hand  side  of 

Eq.  (12)  can  be  approximated  by  a  linear  form  with  an  error  of 
0  (t  -  Tp)2  or  0  (t  -  Tn)2,  respectively.  Then  Eq.  (12)  becomes  a 

quadratic.  When  x  **  x  ,  the  two  roots  have  positive  x^'s;  i.e.,  the 

coning  motion  is  damped.  When  x  **  x  ,  the  solutions  are: 

XR  =  %  (Tn  +  tq0)  +  Js  sgn  (n)  [h  {(m2  +  n2)^  -  m}]'2  (13) 

Tj  =  %  6  ±  %  [h  (  (m2  +  n2)'2  +  m}]  2  (14) 


where 

o  =  (1  -  B)32 

m  =  -  [4  D(xq)  /  (crL)]  +  62  -  (xn  -  T00)  - assumed  >  0 

n  =  2S(x  -  x  )• 
n  oo 

The  lower  signs  in  Eqs.  (13)  and  (14)  yield  negative. Xj,  and  arc  thus 
the  signs  to  use  to  describe  unstable  motion  of  the  gyroscope. 

Eqs.  (13)  and  (14)  are  presented  as  Eqs.  (6-7)  and  (6-8)  in  Ref.  9. 

In  Ref.  9,  Eq.  (6-7)  is  incorrect  in  that  the  factor  "sgn  (n)",  appearing 
in  Eq.  (13),  is  missing  from  the  right-hand  side.  Additionally,  Eq.  (6-6) 
and  the  definition  of  n(=-fr)  are  inconsistent.  Appendix  A  presents  a 
derivation  of  Eqs.  (13)  and  (14). 


E .  Review  of  Assumptions 

It  is  important  to  carry  out  experiments  under  conditions  which  do 
not  violate  the  assumptions  of  the  theory  with  which  measurements  will 
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be  compared.  It  is  now  proper  to  state  or  restate  some  of  the 
important  assumptions  inherent  in  the  theory: 

1.  The  rotor  cavity  is  a  right  circular  cylinder  (no  rounded 

corners) . 

2.  The  unperturbed  motion  of  the  liquid  is  solid-body 
rotation;  also  the  spin  rate  is  large  enough  so  that  in  a  partially- 
filled  cylinder  the  void  forms  a  cylinder  concentric  with  the  rotor 
cavity. 

3.  The  Reynolds  number  of  the  flow.  Re  (=  a is  very 
large.  Viscous  effects  are  thus  confined  to  thin  boundary  layers. 

4.  The  coning  angle  of  the  gyroscope,  |e|,  is  small, 
producing  a  small  perturbation  in  the  solid-body  rotation  of  the  fluid. 

5.  The  gyroscope  motion  is  given  by  Eq.  (6),  and  the 
perturbed  liquid  motion  has  the  same  exponential  dependence  on  time. 

This  implies  that  there  is  an  interval  of  time  in  which  the  motion  of 
the  coupled  gyroscope- liquid  system  is  independent  of  the  realistic 
initial  conditions  but  satisfies  the  small  amplitude  restriction. 

6.  The  coning  frequency  is  in  the  neighborhood  of  a  liquid 
eigenfrequency . 

7.  The  liquid  mass  is  much  smaller  than  the  mass  of  the  gyro¬ 
scope.  This  condition  furnishes  a  simplification  in  the  liquid  moment 
determination. 

8.  The  transverse  moments  of  inertia  about  the  x  -  and  y  - 
axes  are  approximately  equal  to  each  other. 

III.  THE  EXPERIMENTS 


A.  Description 


The  experiments  which  will  be  presented  here  are  those  of  Karpov4, 
Frasier8,  D'Amico10,  Scott  and  D'Amico11,  and  Kitchens  (unpublished). 
All  of  these  experiments  were  performed  at  BRL  using  a  gyroscope  like 
that  described  in  Ref.  4. 

Two  points  to  be  noted  are  that:  (1)  only  one  flexural  pivot  was 
equipped  with  a  strain  gauge,  so  that  only  one  component  of  the  coning 
motion  was  measured;  there  was  no  check  made  of  the  assumption  of 
circular  coning  (see  Section  II-B);  (2)  the  coning  motion  of  the 
empty  gyroscope  was  observed  to  be  damped.  It  was  assumed  that  the 
observed  growth  rate  of  the  coning  motion  of  the  filled  gyroscope 
(Tobs)  could  be  expressed  as  the  sum  of  the  damping  rate  of  the  empty 

gyroscope  (Tempty)  ancl  tlie  growth  rate  of  a  frictionless  filled 
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gyroscope,  i.e. 


T  ,  =  T  .  +  Tt  .  (15) 

obs  empty  I 

This  assumption  appears  reasonable,  although  it  neglects  any  interactive 
effect  between  mechanical  friction  and  presence  of  fluid.  Thus,  the 
experimental  values  of  yaw  growth  rate  recorded  in  this  report  are  x  ^ , 

corresponding  to  the  theoretical  x^  of  Eq.  (14),  not  the  directly 

measured  growth  rates.  Records  indicate  that  a  measured  value  of 

x  ft  =  0.017  s"1  was  used  for  several  of  the  experiments,  so 

empty 

that  x  ^  -  0.324  x  10  4  for  ft  =  523.6  rad/sec,  the  value  used  in  the 

empty 

experiments . 

The  electrical  signal  from  the  strain  gauge  was  recorded  on  an 
oscillograph.  Since  the  output  of  the  strain  gauge  was  linear  with 
angular  displacement  of  the  flexural  pivot,  the  oscillograph  record 
yielded  a  plot  of  coning  angle  versus  time.  A  set  of  amplitude  readings 
from  the  oscillograph  record  was  plotted  logarithmically  against  time. 
These  points  fell  about  a  straight  line  (Eq.  (6a)),  until  the  amplitude 
reached,  typically,  one  degree.  The  observed  yaw  growth  rate  was  deter¬ 
mined  as  the  slope  of  this  line.  The  frequency  x^  of  the  coning  motion 

was  measured  by  counting  the  number  of  oscillations  between  two  known 
times  on  the  oscillograph  record.  The  rotational  speed  of  the  cylinder 
was  measured  with  a  stroboscope  * 

At  this  late  date  it  is  difficult  to  ascertain  the  errors  of  past 

measurements.  A  study  of  these  experiments  yields  estimated  relative 

errors  of  2%  and  3%  for  x  and  x  ,  respectively.  These  are  the  errors 

R  1 

indicated  by  the  error  bars  in  the  figures  in  this  section;  the  apparent 
scatter  in  the  data  seems  to  be  consistent  with  these  estimates. 
Estimated  errors  for  cylinder  dimensions  and  fluid  parameters  are  as 
follows : 


Aspect  Ratio  of  Cylinder  (c/a)  %-l% 

b2/a2  (or  d2/a2)  2% 

Reynolds  Number  (Re)  5% 

Liquid  Density/ Axial  Moment  of  Inertia  (p/L)  1% 


Except  for  the  error  in  aspect  ratio,  the  effects  of  these  errors  on 
computation  of  x^  are  small.  The  effect  of  the  error  in  c/a  can  be 

significant,  as  discussion  in  Section  III-B  will  demonstrate. 
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There  is  a  division  between  Karpov's  experiment  and  the  others,  in 
that  Karpov  used  constant  values  of  L  and  x  ,  varying  the  fill  ratio  of 

the  cylinder  (hence  the  liquid  eigenfrequency) ,  while  the  other  experi¬ 
ments  varied  L  and  x^  while  maintaining  constant  fill  ratio.  Karpov 

made  separate  measurements  of  x  and  xT  as  functions  of  b2/a2.  The 

K  1 

others  made  simultaneous  measurements  of  td  and  xT  as  functions  of  L  and 

K  1 

x^,  though  unfortunately  they  did  not  measure  or  record  L  or  This 

latter  fact  prevents  a  full  comparison  of  their  experimental  data  with 

the  theory.  The  nature  of  the  comparison  which  can  be  made  is  indicated 
in  Section  III-C.  The  Reynolds  number  was  varied  by  using  liquids  of 
varying  viscosity,  keeping  the  cylinder  speed  fixed  at  523.6  rad/s.  This 
means  that  there  was  no  check  of  the  effect  of  gravity  in  the  experi¬ 
ments  with  the  partially-filled  cylinder.  Details  of  cylinder  dimen¬ 
sions,  liquid  properties,  moments  of  inertia,  etc.,  and  the  measured 

values  of  x  and  xT  are  tabulated  in  Appendix  B. 

K  1 

B.  The  Results  of  Karpov 

Figs.  3  and  4  present  x  and  -xT  for  the  mode  n  =  j  =  1  plotted 

K  1 

against  b2/a2  for  the  two  liquids  that  Karpov  used,  labeled  #1  (Re  = 

5.2  x  105)  and  #2  (Re  =  5.2  x  103).  The  experimental  values  are  indicated 
by  circles  and  diamonds.  The  theoretical  values,  solid  curves,  are 
obtained  from  Eqs.  (13)  and  (14).  The  plots  of  -t  against  b2/a2  are 

analogous  to  Figs.  4  and  5  in  Ref.  4,  Fig.  1  in  Ref.  6,  and  Fig.  6-2  in 
Ref.  9,  where  the  abscissa  is  the  parameter  x  .  The  significant  differ¬ 
ence  between  these  figures  and  Figs.  3  and  4  is  that  the  latter  two 
contain  true  theoretical  curves,  as  opposed  to  semi-empirical  curves 
in  the  former  (discussed  in  Section  I). 

The  predicted  trend  of  xD  seems  to  appear  in  the  data,  as  does  the 

K 

"jump"  at  b  "/a2  «  0.15  for  liquid  #1,  However,  the  data  are  consistently 
about  3%  higher  than  the  theoretical  predictions. 

For  the  Re  =  5.2  *  105  case,  the  measured  values  of  x ^  agree  with 

the  predicted  values  to  within  the  estimated  error.  For  the  Re  = 

5.2  x  103  case,  the  measured  peak  value  of  x^  is  about  10%  greater  than 

the  prediction,  and  there  is  an  apparent  shift  of  the  measured  curve  to 
the  left  of  the  theoretical  curve.  A  study  of  the  sensitivity  of  x ^  to 

small  changes  in  aspect  ratio  was  made  with  the  intention  of  determining 
whether  this  shift  could  be  explained  by  uncertainties  in  aspect  ratio 
measurements.  In  Figure  5,  x ^  vs  b2/a2  curves  are  shown  for  c/a  greater 

and  less  than  the  experimental  aspect  ratio  by  1/2  percent.  For  the 
less  viscous  fluid,  the  curves  bracket  the  data;  the  apparent  discrep¬ 
ancy  between  the  experimental  and  theoretical  peaks  is  not  significant. 
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For  liquid  #2,  the  conclusion  is  not  as  definite.  A  positive  correction 
of  1/2  percent  to  the  aspect  ratio  brings  the  theoretical  and  experi¬ 
mental  values  of  Xj  into  closer  agreement,  including  the  positions  of 

the  peaks;  however,  the  theoretical  curve  exhibits  considerably  more 
asymmetry  about  the  peak  than  do  the  data. 

C .  Results  of  D'Amico,  Kitchens,  and  Frasier 

The  quantities  t^  and  L  were  not  recorded  for  these  experiments. 

Fortunately,  the  ring  combinations  used  were  recorded,  and  all  of  the 
rings,  rotor,  and  cylinder  are  extant,  so  that  measurements  could  be 
made  of  L.  However,  T  (or  equivalently,  =  L/T)  corresponding  to  the 

experimental  conditions  could  not  be  determined  because  the  location  of 
the  counterweight  used  to  adjust  the  position  of  the  center  of  gravity 
of  the  gyroscope  was  not  recorded. 

In  a  particular  set  of  experiments,  L  and  were  varied  for  each 

case,  so  that  t  =  t  (L) .  Since  the  t  's  in  these  experiments  were  not 
n  n  n 

recorded,  tD  and  x,  could  not  be  evaluated  by  Eqs.  (13)  and  (14) 

explicitly  in  terms  of  measured  t  .  The  relationship  x  (L)  was  produced 

indirectly  by  solving  Eq.  (13)  numerically  for  rn  for  all  pairs  of 

measured  L  and  xR,  and  fairing  curves  through  the  plotted  xn  vs  L 

points.  In  most  cases,  one  straight  line  fitted  the  points  quite  well, 
the  r.m.s.  error  of  L  about  the  line  being  about  1/2%  of  the  mean. 

These  points  correspond  to  changes  in  the  increment  rings  only  (see 
Appendix  B) .  In  two  experiments,  however,  the  main  ring  was  changed  to 
extend  the  range  of  L,  producing  a  discontinuity  in  the  x^  vs  L  curve; 

in  these  cases  the  points  were  fitted  with  two  straight  lines,  one  for 

each  main  ring.  This  fit  led  to  discontinuities  in  the  curves  of  x  vs 

x  (Figs .  6  and  7) . 

K 

Values  of  x„  and  x  were  then  evaluated  by  Eqs.  (13)  and  (14)  for 
R  I 

each  series  of  L's  employed  in  the  experiments,  and  plots  of  Xj  vs  tr 

were  drawn,  which  could  be  compared  directly  with  experiment.  These 
plots  do  not  provide  a  total  check  of  the  theory,  since  part  of  the 
input  (t  )  was  based  on  the  prior  assumption  of  the  theory's  validity. 

It  was  stated  in  Section  I  that  Frasier's  comparison  of  his  data8 

to  the  theory  was  flawed;  the  nature  of  the  flaw  is  now  clear.  Since 

neither  r  nor  L  was  measured  during  his  experiments  (except  that  a 
n 

single  "nominal"  L  was  noted),  he  could  not  have  correctly  evaluated 
Eq.  (14)  to  find  t , .  It  seems  likely  that  he  did  what  later  experi¬ 
menters  did;  namely,  use  xn  in  place  of  x  in  Eq.  (14)  . 

K  n 
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Figs.  6  to  10  present  plots  of  -tt  vs  td  for  all  of  the  experiments, 

1  R 

including  four  different  liquids  at  79%  fill,  two  liquids  at  100%  fill, 

and  three  liquids  in  a  rodded  cylinder.  The  error  bars  indicate  ±  2% 

error  in  t  and  ±  3%  error  in  :  ,  General  comments  about  the  agreement 
K  1 

between  theory  and  experiment  are  that:  (1)  the  amplitude  at  the  peak 
value  of  Tj  agrees  to  within  10%  when  Re  >  4  x  104;  (2)  the  amplitude 

at  the  peak  agrees  within  25%  at  lower  Reynolds  numbers;  (5)  in  general, 
the  theory  predicts  a  more  rapid  decay  away  from  the  peak  than  was 
observed,  particularly  to  the  right  of  the  peak.  In  most  of  the 
figures,  there  is  a  suggestion  of  a  shifting  of  the  predicted  peak 
response  to  the  left  or  to  the  right  of  the  observed  peak.  However, 

Fig.  11,  showing  the  effect  of  a  ±  1/2%  change  in  aspect  ratio  on  two 
of  the  curves,  demonstrates  that  the  apparent  shifts  are  not  significant. 


IV.  CONCLUDING  REMARKS 

Renewed  emphasis  is  currently  being  placed  on  research  into  the 
dynamics  of  liquid-filled  shell,  in  particular  by  means  of  gyroscope 
experiments.  In  the  planning  of  future  work,  it  is  necessary  to  have  a 
clear  picture  of  what  has  been  accomplished  in  the  past.  Examination 
of  available  records  and  reports  has  revealed  inadequacies  in  the 
results  of  previous  investigators,  such  as  incomplete  recording  of  data 
(e.g.,  moments  of  inertia)  and  misleading  comparisons  of  theoretical 
and  experimental  results.  This  report  was  written  in  an  effort  to 
provide  a  correct  and  unified  picture  of  past  work. 

Comparisons  show  that  predicted  and  measured  frequencies  of 
greatest  instability  of  motion  agree  to  within  the  experimental  error; 
the  growth  rates  at  those  frequencies  agree  to  within  10%  for  Re  >  104, 
and  to  within  25%  for  lower  Reynolds  number.  All  work  discussed  here 
was  performed  with  the  liquid  completely  spun-up;  i.e.,  in  rigid  body 
rotation.  The  conclusions  in  Ref.  9  based  on  this  work  are  valid;  namely, 
that  free-oscillation  frequencies  of  the  liquid  are  well-predicted  by 
the  Stewartson-Wedemeyer  theory,  and  that  coning  frequencies  and  yaw 
growth  rates  near  resonance  agree  qualitatively  with  theoretical 
predictions.  However,  the  experiments  on  which  these  conclusions  are 
based  cover  only  a  small  range  of  coning  frequencies,  ratios  of  liquid 
mass  to  solid  mass,  and  aspect  ratios,  and,  additionally,  only  a  single 
rotation  speed.  Also,  no  direct  measurement  has  been  made  of  the 
motion  of  the  liquid.  Therefore,  further  work  is  required  in  this  area 
to  obtain  a  more  distinct  quantitative  picture  of  the  effect  of  the 
liquid  on  the  stability  of  the  shell. 


20 


OUTER 
GIMBAL  PIVOT 


INNER 

/■*  I  A  A  n  A  t 


ROTOR 

CAVITY 


GIMBAL  PIVOT 


Figure  1.  Schematic  Diagram  of  Gyroscope 


Figure  2.  Schematic  Diagram  of  Coordinates 
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Figure  3.  Coning  Frequencies  and  Yaw  Growth  Rates  for  Karpov  Liquid 

#1  Experiments  (a  =  3.150  cm,  c  =  9.691  cm,  p  =  0.818  g/cm3, 
L  =  1.968  x  105  g  cm2,  v  =  0.01  cm2/s,  tn  =  0.0534; 

n  =  1 ,  j  =  1 ) 
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Figure  4.  Coning  Frequencies  and  Yaw  Growth  Rates  for  Karpov  Liquid 

#2  Experiments  (a  =  3.150  cm,  c  =  9.691  cm,  p  =  0.818  g/cm3, 
L  =  1.968  x  105  g  cm2,  v  =  1.00  cm2/s,  T  =  0.0534; 

n  =  1,  j  =  1) 
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Figure  5.  Sensitivity  of  Yaw  Growth  Rate  to  Change  in  Cylinder  Aspect 
Ratio  in  Karpov's  Experiments  (a  =  3.150  cm,  c  =  9.691  cm, 

L  =  1.968  x  105  g  cm2,  t  =  0.0534;  n  =  1,  j  =  1) 
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Figure  6.  -  Tj  vs  for  D'Amico's  79%  Filled  Cylinder 

Experiments,  Re  =  5.2  x  105  and  Re  =  4.0  x  104 
(a  =  3.153  cm,  c  =  9.500  cm,  (b/a)2  =  0.210; 
n  =  1 ,  j  =  1) 
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Figure  7.  -  t.  vs  xR  for  D'Amico's  79%  Filled  Cylinder 

Experiments,  Re  =  1.1  x  104  and  Re  =  5.2  x  103 
(a  =  3.153  cm,  c  =  9.500  cm,  (b/a)2  =  0.210; 
n  =  1,  j  =  1) 
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Figure  8.  -  xT  vs  xD  for  100%  Filled  Cylinders,  Experiments  of 

I  K 

D'Amico  and  Kitchens  (a  =  3.153  cm,  c  =  9.928  cm; 
n  =  1 ,  j  =  1 ) 
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Figure  9.  -  t.  vs  tr  for  Frasier  (Rod)  Experiments,  Re  =  5.2  x  105 

and  Re  =  1.70  x  105  (a  =  3.153  cm,  c  =  9.030  cm, 

(d/a)2  =  0.023;  n  =  1,  j  =  1) 
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Figure  10. 


-  t T  vs  td  for  Frasier  (Rod)  Experiments,  Re  =  4.0  *  104 


(a  =  3.153  cm,  c  =  9.030  cm,  (d/a)2  =  0.023;  n  =  1, 
j  =  1) 
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(a) 


(b) 


Figure  11.  Sensitivity  of  -  ij  vs  Curves  to  Change  in  Cylinder 

Aspect  Ratio:  (a)  D'Amico  100%  Filled  Cylinder, 
a  =  3.153  cm,  c  =  9.928  cm,  j  =  1,  n  =  1; 

(b)  Frasier  Cylinder  With  Inner  Burster,  a  =  3.153  cm, 
c  =  9.030  cm,  (d/a)2  =  0.02,  n  =  1,  j  =  1 
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LIST  OF  SYMBOLS 


a 

b 

c 

d 

D 

j 

L 

M 

o 

n 

R 

Re 

t 

T 

x,  y,  z 

x/  >  y' , 

a 

3 

6 

6a,  5c, 6d 

r 

* 

r 


cross-sectional  radius  of  cylindrical  cavity  (cm) 

radial  coordinate  of  free  surface  in  partially-filled 
cylinder  (cm) 

half-height  of  cylindrical  cavity  (cm) 
cross-sectional  radius  of  inner  burster  (cm) 
liquid  moment  residue  function,  Eq.  (7)  (g  cm  ) 
index  of  axial  perturbation  mode 

axial  moment  of  inertia  of  empty  gyroscope  (g  cm  ) 

moment  function  due  to  displacement  of  center  of  mass  from 
pivot  point,  Eq.  (1)  (g  cni  /s2) 

index  of  radial  perturbation  mode 

liquid  moment  residue  function,  Eq.  (8)  (non-dimensional) 

(=  afi2/v)  Reynolds  number  (non-dimensional) 
time  (s) 

transverse  moment  of  inertia  of  empty  gyroscope  about  axis 
through  pivot  point  (g  cm2) 

laboratory  fixed  rectangular  coordinates,  z-axis  coinciding 
with  unyawed  cylinder  axis.  Figs.  1  and  2  (cm) 

z'  non-rotating  rectangular  coordinates,  z' -axis  coinciding 
with  yawing  cylinder  axis.  Figs.  1  and  2  (cm) 

coning  angle.  Fig.  2  (radians) 

=  4  MQT/(Lfi)2  (non-dimensional) 

viscous  perturbation  decay  rate/Sl  (non-dimensional) 

effective  changes  to  a,  c,  and  d  due  to  boundary  layers, 

Eq.  (9)  (cm) 

liquid  moment  function  of  x  in  Eq.  (3)  (g  cm2/s2) 

2  2 

liquid  moment  function  in  Eq.  (1)  (g  cm  /s  ) 
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LIST  OF  SYMBOLS  (continued) 


e(=  e  +ie  ) 

x  y 

v 

P 

a 


complex  yaw  (0  and  0^  being  x  and  y  components, 
x  y 

respectively)  (non-dimensional) 

2 

kinematic  viscosity  of  liquid  (cm  /s) 

3 

density  of  liquid  (g/cm  ) 

5  ci  - 


T(=  TR+iTp 

T 

O 


oo 


ov 


T 

P 


exponent  term  in  exponential  solution,  Eq.  (2) 
inviscid  perturbation  eigenfrequency/ft  (non-dimensional) 
viscous  perturbation  eigenfrequency/ ft  (non-dimensional) 

“  Too  +  1  Eq* 

-  (yaw  growth  rate)/Sl  (non-dimensional) 

nutational  frequency  of  empty  gyroscope/Q  ,  Eq.  (4) 
(non-dimensional) 

precessional  frequency  of  empty  gyroscope/ft  ,  Eq.  (4) 
(non-dimensional ) 

coning  frequency  of  liquid-filled  gyroscope  (non- 
dimensional) 


SI 


spin  rate  of  cylinder  (radians/s) 
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APPENDIX  A:  DERIVATION  OF  EQUATIONS  (13)  AND  (14) 


We  start  with  Eq.  (6-2)  of  Ref.  9: 


T  =Cl/2)(Xn  +  Tov)  ±  (  Rl/2)  (Tn  -  Tov)  ]2  -  [  D|  /  (CTL)  > 


1/2 


Recalling  that  =  xqq  +  i5,  Eq.  (A.l)  can  be  written 
t  =(.1/2)  (Tn  +  XQo)  +(l/2) i<5  ±(l/2)i  ( (4  |  D | )  /  (crL)  +  6 2  -  (xn~  x 


+  2i6  (t  -  t  )  ) 
n  00 


1/2 


If  we  now  define 

m  =  (4 1  D | )  /  (aL)  +  62  -  (xR  -  tqo)2 

Or 

and  n  =  2  6  (x  -  x  )  > 

n  00 

where  m>0  for  the  ranges  of  parameters  considered.  Eq.  (A. 2) 
x  =U/2')  (xn  +  xqo)  +(^l/2)i6  ±(1/2  i  (m  +  in)1^2 

Or 

This  corresponds  to  Eq.  (6-6)  in  Ref.  9,  with  n  replacing  n. 
Eq.  (6-6)  of  the  reference  is  in  error,  given  the  definition 
there . 

.  ^  i0 

Let  m  +  1  n  h  X  e  , 

where 

Or  Or 

A2  =  m2  +  n2  ,  6  =  tan-1  (n/m)  ,  sgn(0)  =  sgn(n)  ,  -  if  <  0 

Then 

±  i  (m  +  i  n)^2  =  ±  A^2  [-  sin  (0/2)  +  i  cos  (0/2)] 
For  -  tt/2  <  0/2  <  tt/2  , 

sin  (0/2)  =  sgn(0)  [(1  -  cos  0)/2]^2  , 

and  cos  (0/2)  =  [(1  +  cos  0)/2]^2 


(A.l) 


(A. 2) 


can  be  written 


(A. 3) 


Note  that 
if  n 


(A. 4) 


<  it  .  (A. 5) 
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Thus 


±  i  (m  +  i  n)^2  =  +  x^2  [-sgn(0)  {(1  -  cos  0)/2}^2 

+  i  {(1  +  cos  6)/2}1/2]  .  (A. 6) 

Using  the  following  trigometric  identity: 

cos  (tan-1  x)  =  (1  +  x2)-^2  for  -  ir/2  <  tan-1  x  <  tt/2  , 
and  Eq.  (A. 5),  we  get 

±  i  (m  ♦  i  S’)1'2  -  ±  2-1''2  [-sg„(Sj  {(>2  ♦  -  m}1'2 

*  i  { (m2  *  n2)1/2  +  m)1/2]  .  (A. 7) 

Inserting  Eq.  (A. 7)  into  Eq.  (A. 3)  and  separating  real  and 
imaginary  parts,  we  get 

xR  =vl/2)(rn  +  tqo)  +  (l/2) sgn(n)  [{  (m2  +  n2)1^2  -  m}/2]1//2 

Tj  =(1/2)6  ±(l/2)  [{(m2  +  n2)1/2  +  m}/2]1/2  ,  (A. 8) 

which  are  Eqs.  (13)  and  (14)  in  this  report. 
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APPENDIX  B.  DATA  FROM  EXPERIMENTS 


For  the  convenience  of  future  investigators,  relevant  data  from 
the  experiments  described  in  this  report  are  tabulated  here.  All 
frequencies  and  growth  rates  listed  here  were  obtained  by  the  original 
researchers  from  the  raw  data;  no  re-reduction  of  raw  data  was 
performed  by  the  present  authors,  except  to  verify  the  reduction  methods. 

In  all  experiments  except  those  of  Karpov,  the  axial  moment  of 
inertia  of  the  gyroscope  was  varied  by  changing  rings  mounted  on  the 
rotor.  These  rings  were  of  two  types:  main  rings,  in  five  sizes;  and 
increment  rings,  in  four  sizes.  At  any  one  time,  only  one  main  ring 
was  used,  but  any  or  all  increment  rings  could  be  used.  The  increment 
rings  were  mounted  above  or  below  the  main  rings,  the  two  positions 
being  distinguished  as  "top"  and  "bottom”  in  the  experimental  records. 

The  axial  moments  of  inertia  of  the  components  of  the  gyroscope  were 
measured  by  the  present  authors  and  are  shown  in  Table  B-l,  together 
with  the  corresponding  masses.  There  were  three  increment  rings  of 
size  .001  and  two  of  each  of  the  other  sizes.  The  values  shown  in 
Table  B-l  represent  mean  values  of  the  rings  of  each  size. 


TABLE  B-l.  AXIAL  MOMENT  OF 
Component  Mass 


INERTIA  OF  GYROSCOPE  COMPONENTS 

(g)  Axial  Moment  of 

Inertia  (g  cm^) 


Steel  Rotor 

2554 

3.221 

Plastic  Cylinder 

933 

0.841 

Main  Ring 

i 

1650 

6.238 

Main  Ring 

2 

3160 

13.87 

Main  Ring 

3 

4455 

21 . 74 

Main  Ring 

4 

6510 

31.03 

Main  Ring 

5 

7725 

39.22 

Increment 

Ring  . 

001 

180 

0.48 

Increment 

Ring  . 

002 

309 

0.93 

Increment 

Ring  . 

0025 

367 

1.10 

Increment 

Ring  . 

005 

615 

2.12 

In  the  experimental  records,  numbers  are  given  which  indicate  the 
total  ’’value"  of  the  increment  rings  used.  For  example,  in  one  case 
the  increments  might  be  recorded  as  ’’top-.  003,  bottom- .  005” .  This 
could  represent  any  of  the  combinations  shown  in  Table  B-2.  Unfor¬ 
tunately,  not  all  of  the  combinations  have  the  same  axial  moment  of 
inertia.  In  this  work,  wherever  there  were  several  possible  increment 
combinations,  the  mean  moment  of  inertia  of  the  combinations  was  used. 
In  the  worst  case,  this  led  to  an  inaccuracy  of  less  than  one  percent 
of  the  total  axial  moment  of  inertia  of  the  gyroscope. 
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TABLE  B-2.  INCREMENT  RING  COMBINATIONS  WHICH  COULD  BE  LISTED  AS 

"TOP- . 003,  BOTTOM-. 005M 

Top  Increments  Bottom  Increments  Total  Axial  M.O.I.  ^ 

of  Increments  (g  cm") 


.001  , 

.002 

.005 

3.53 

X 

10 

.001  , 

.002 

2  x  .0025 

3.62 

X 

10 

3  x 

.001 

.005 

3.56 

X 

10 

3  x 

.001 

2  x  .0025 

3.65 

X 

10 

Mean 


3.59  x  10 


4 


In  the  following  tables,  the  data  describing  each  experiment  are 
presented.  The  date  indicated  for  each  experiment  is  approximate. 


TABLE  B-3 .  THF,  EXPERIMENT  OF  KARPOV  (1965) 


a  =  3.150  cm 

t  =  0.0534 

n  2 
v  =  0.01  cm  /s~ 

p  =  0.818  g/cnf 


c  =  9.691  m  _ 

L  =  1.968  x  10'  g  cm 
2 

v  =  1.0  cm  /s-7 
p  =  0.968  g/cm' 


,2  .  2 
b  /a 

tr 

"TI 

tr 

-TI 

0 

.0540 

A 

.0552 

.439 

.0248 

.0542 

.038  x  10 

.0552 

.  575 

.0496 

.  0543 

.076 

.0548 

.707 

.0744 

.0545 

.  115 

.0549 

.993 

.0992 

.0545 

.153 

.0549 

1 . 394 

.1240 

.0544 

.458 

.0552 

2.024 

.  1322 

.0545 

.917 

.0554 

2.082 

.  1405 

.0544 

1.872 

.0551 

1.967 

.1446 

3.400 

.1488 

7.563 

.0552 

1.872 

.1504 

9.015 

.1537 

9.855 

.1554 

9.969 

.1570 

.0552 

5T- 

.1587 

8.747  » 

l 

! 
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2  2 
bZ/a 

TABLE 

i 

! 

B-3. 

tr 

THE  EXPERIMENT 

"TI 

OF  KARPOV  (1965) 

tr 

(continued) 

-Ti 

.1620 

I 

7.219 

x  10"4 

.1653 

5.691 

. 

.1686 

2.139 

.1736 

.0551 

1.070 

.1818 

.0550 

.477 

.0555 

.859  x  10 

.1983 

.0552 

.497 

.2066 

.0548 

.095 

.2231 

.0552 

.420 

.2314 

.0548 

.038 

T 

TABLE 

B-4.  THE  EXPERIMENT  OF  D'AMICO 

(1969) 

a  =  3. 

153  cm 

c  =  9.500 

cm  b^/a 

2 

=  0.21 

V  = 

2 

0.01  cm  /s  , 

p  =  0.818 

g/cm3 

Main 

Increment 

Rings 

2 

tr 

-TI 

Ring 

Top 

Bottom 

L(g  cm  ) 

0 . 0  0  0  0 

0.0010 

1 . 878E  +  05 

0. 0498 

1 . 222E-04 

0. 0010 

0. 0010 

1 . 926E+05 

0. 0507 

1 . 375E-04 

-| 

0 . 0020 

0. 0010 

1 . 973E+05 

0.0515 

1 . 623E-04 

■“i 

0  .  0  0  2  0 

0. 0020 

2.01 8E+05 

0 . 0522 

1 . 891E-04 

2 

0. 0025 

0. 0025 

2. 051E+05 

0.0532 

4 . 2  5  9  E  ~  O  4 

p 

0 .  0  0  2  5 

0 . 0035 

. 099E+05 

0. 0542 

5 . 806E-04 

-i 

0 . 0035 

0 . 0035 

2. 147E+05 

0. 0549 

7 . 6  7  8  E  ~  0  4 

2 

0. 0045 

0 .  0  0  3  5 

2. 1 93E+05 

0. 0558 

8. 537E-U4 

2 

0 . 0040 

0 . 0050 

2. 234E+05 

0. 0565 

y .  3  4  6  E — 0  4 

0 . 0050 

0 . 0050 

. 276E+05 

0. 0574 

6 . 2  2  6  E  —  U  4 

0 . 0050 

0 . 0060 

. 320E+05 

0. 0582 

4 . 6  9  8  E  ~  U  4 

0 . 00b 0 

0 . 0  0  6  0 

2. 36SE+05 

0. 0589 

4 . 736E-04 

-• 

0 . 0  0  6  0 

0 . 0070 

.  4  U  2  E  +  0  5 

0 . 0600 

2. 327E-Q4 

•-i 

0 . 0070 

0 . 0  0  7  0 

2. 443E+05 

0. 0610 

2. 215E-04 

0 . 0  0  7  5 

0 . 0  075 

. 494E+05 

0. 0614 

1 . 719E-04 

0 . 0085 

0 . 0  0  7  5 

2. 537E+05 

0. 0623 

1 . 700E-04 

0 . 0  0  8  5 

0 . 0095 

.  6 1 7  E + W  5 

0 . 0636 

1 . 490E-04 

0 . 0  1  0  5 

0. 0095 

2. 71 1E+05 

0. 0651 

1 . 889E-04 

39 


TABLE  B-4 .  THE  EXPERIMENT  OF  D'AMICO  (1969)  (continued) 


2 

v=0.13cra/s  p  = 

0.940  g/cm3 

Main 

Increment  Rings 

2 

T 

-  T 

Ring 

Top 

Bottom 

L (g  cm  ) 

R 

I 

0 . 8000 

0 . 0020 

1 . 925E+05 

9. 0460 

8. 594E-05 

ki 

0 . 0020 

0.0020 

2.01 8E+05 

0 . 043 U 

1 . 089E-04 

2 

0 . 0060 

0 . 0000 

2. 100E+05 

0. 0490 

1 . 509E-U4 

2. 

0 . 0060 

0 . 0020 

2. 189E+05 

0. 0510 

2. 177E-04 

2 

0 . 0050 

0 . 0050 

2 . 2  7  6  E  +  0  5 

0. 0530 

2. 597E-04 

2 

0 . 0070 

0 . 0050 

2. 365E+05 

0.0540 

3 . 686E-04 

u— 

0 . 0070 

0 . 0070 

2. 448E+05 

0 . 0550 

3. 992E-04 

lL 

0 . 0075 

0 . Q035 

2. 537E+05 

0. 0560 

4 . 030E-04 

,3 

0 . 0095 

0 . 0085 

2. 617E+05 

0 . 0580 

3. 304E-04 

£l 

0 . 0  1  0  5 

0 . 0095 

2. 71 1E+05 

0 . 0590 

2. 941E-04 

3 

0 . 0000 

0 . 0000 

2 . 620E+05 

0. 0620 

1 . 776E-U4 

3 

0 . 0000 

0. 0020 

2. 715E+05 

0 . 0630 

1 . 261E-04 

3 

0 . 0020 

0 . 0020 

ii .  8  U  0  E  t  U  5 

0. 0650 

9. 931E-05 

■-* 

0 . 0000 

0 . 0060 

2. S90E+05 

0.067O 

8. 021E-05 

v  =  0.49 

2, 

cm  /s  ,  p  = 

0.960  g/cm3 

Main 

Increment  Rings 

Ring 

Top 

Bottom 

Ug  cm  ) 

tr 

TI 

0. 0020 

0 . 0020 

2 . 0 1  8  E  +  0  5 

0. 0460 

7. 257E-05 

0. 0010 

0 . 0050 

2. 100E+05 

0. 0480 

8 .  i'  8  5  E  0  w1 

/  1 

0. 0030 

0 . 0050 

2. 189E+05 

0 . 0490 

1 . 1 65E-04 

0. 0050 

0 . 0050 

2 . 2 7 6 E + U  5 

G. 0510 

1 . 3 37 E- 04 

0. 0070 

0 . 0050 

2 . 3 6  5  E + 0  5 

0 . 052 G 

1 . 604E-G4 

0. 0070 

0. 0070 

2 . 4  4  8  E + 0  5 

0 . 0540 

1 . 929E-04 

0. 0085 

0. 0075 

2. 537E+05 

G . 0550 

2 . 06  3 E - 0  4 

0. 0085 

0. 0095 

2. 61 7E+05 

0 . 0560 

2. 254E-04 

0 . 0 1 05 

0. 0095 

2.71 1E+05 

0 . 0  5  8  0 

2. 655E-04 

0. 0000 

0 . 0  0  0  0 

2 . 620E+05 

0 . 0  6  0  0 

2. 483E-G4 

J 

0. 0000 

0 . 0020 

2. 715E+05 

0. 0620 

2. 807E-04 

3 

0 . 0020 

0 . 0020 

2 . 8  0  3  E  +  0  5 

0. 0630 

2. 426E-04 

0 . 0000 

0. 0035 

2. 778E+05 

0. 0630 

2 . 082 E- 04 

~l 

0 . 0020 

0 . 0035 

2. 873E+05 

0 . 0640 

2. 024E-04 

0 . 0  0  4  0 

0 . 0035 

2. 966E+05 

0. 0660 

1 . 7  7  6  E  -  0  4 

•  J 

0 . 0030 

0 . 0070 

.  U  7  3  E  +  0  5 

0 . 0680 

1 . 509E-G4 

0 . 0070 

0 . 007 0 

3. 238E+05 

0. 0720 

1 . 337 E- 04 

•~t 

0.0095  • 

0 . 0105 

3. 5O1E+05 

0 . 0750 

9. 358E-05 
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TABLE  B-4 .  THE  EXPERIMENT  OF  D’AMICO  (1969)  (continued) 

2  z 

v  =  1.0  cm  /s  ,  p  =  0.968  g/cm' 


Main  Increment  Rings 

R j ng  Top  Bottom _  L(g  cm2)  _ ^R _ ~TI 


0 . 0000 

0. 0000 

1 . S30E+05 

0. 0430 

4. 5S4E-05 

2 

0. 0020 

0. 0020 

2. 01SE+05 

0.0460 

3. 212E-05 

2 

0. 0045 

0 . 0  0  3  5 

2. 193E+05 

0. 0490 

9. 549E-05 

2 

0 . 0060 

0 . 0060 

2 . 363E+05 

0. 0520 

1  .  1  3  4  E  -  0  4 

2 

0 . 0  0  7  5 

0 . 0035 

2. 537E+05 

0 . 0530 

1 . 5 35 E- 04 

2 

0 . 0095 

0. 0105 

2. 71 1E+05 

0. 0610 

1 . 719E-04 

o 

0. 0000 

0 . 0  0  2  0 

2. 715E+05 

0. 0620 

1 . 391E-U4 

o 

0. 0020 

0. 0020 

2. 303E+05 

0. 0640 

1 . 3 1  4  E  -  0  4 

•-» 

0. 0025 

0 . 0035 

2. SS9E+05 

0. 0650 

1 . 7  7  6  E — 0  4 

2 

0. 0045 

0 . 0035 

2. 9S3E+05 

0. 0660 

1 . 642E-04 

0 . 0050 

0 . 0050 

3. 066E+05 

0. 0630 

1 . 623E-04 

o 

N-' 

Q . 0050 

0. 0070 

3. 155E+05 

0. 0690 

1 . 413E-04 

V1 

0. 0070 

0. 0070 

3. 23SE+05 

0. 0710 

1 . 3  3 1  E  U  4 

3 

0. 0105 

0. 0095 

3. 501 E+05 

0.0750 

1 . 070E-04 

TABLE 

B-5.  THE 

EXPERIMENT  OF 

D’AMICO  (1974) 

a  =  3. 

.153  cm 

c 

=  9.928  cm 

2  2 
b  /  k  = 

0 

Main 

v  = 

0.01  cm2/ s 

,  p  =  0 

.818  g/ci2 

Ring 

Top 

Bottom 

L(g  cm2) 

tr 

"TI 

0 .000 0 

0 . 0000 

1 . S30E+05 

0.0450 

5.  7  3  U  E  -  0  4 

2 

0. 0010 

0. 01000 

1 . 378E+05 

0. 0456 

7  .  0  6  6  E  —  0  4 

cL 

0. 0010 

0. 0010 

1 . 926E+05 

0 . 047  6 

9. 773E-04 

2 

0 . 002  0 

0. 0010 

1 . 9  7  3  E  +  05 

0.0434 

3 . 3  4  6  E  -  0  4 

2 

0 . 0020 

0 . 0020 

2.01 SE+05 

0.0493 

3. 862E-04 

o 

0 . 0020 

0 . 0030 

2 . 0  6  6  E  +  0  5 

0.0504 

7. 31  IE-04 

2 

0 . 0  0  2  0 

0 . 0  01  0!  0 

2 .  Li  6  6  E  +  0i  5 

0 . 0506 

7 . 023E-04 

2 

0 . 0  0  3  0 

0 . 0030 

2. 1 12E+05 

0.0514 

5 . 500E-04 

0. 0035 

0. 0035 

2. 147E+05 

0. 0523 

3. 304E-04 

2 

0.0035 

0. 0045 

2. 193E+05 

0.0523 

3. 151E-04 

2 

0 . 0  0  5  - 

0 . 0050 

2 . 276E+05 

0.0533 

2. 292E-04 

2 

0 . 0  0  5  0 

0 . 00  6  01 

2 . 320E-+05 

0. 0553 

1 . S14E-04 

2 

0 . 006  0 

0 . 01060 

2. 36SE+05 

U . 0564 

1 . 757E-04 

2 

0 . 0  01  6  01 

0 . 0  0  7  0 

2 . 402E+05 

0 . 0573 

1 . 3  7  5  E  -  0  4 

2 

0 . 0070 

01 . 0070 

2. 443E+05 

0. 0533 

1 . 039E-04 

2 

0 . 0  0  7  01 

0 . 0030 

2. 496E+05 

0.0593 

1 . 31SE-04 
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TABLE 

B-6.  THE 

EXPERIMENT  OF 

KITCHENS  (1976) 

a  =  3. 

153  cm 

c 

=  9.928  cm 

.2.  2 
b  /a  = 

0 

v  = 

0.58  cm2/s 

,  P  = 

0.960  g/cm'^ 

Main 

Ring 

Increment  Rings 

Top  Bottom 

L(g  cm2) 

tr 

"TI 

2 

0.0000 

0.0010 

1 . 878E+05 

0.0437 

1.53E-04 

2 

0.0000 

0.0020 

1 . 925E+05 

0.0447 

1.39E-04 

2 

0.0010 

0.0020 

1 . 973E+05 

0.0456 

1.45E-04 

2 

0.0020 

0.0020 

2 . 018E+05 

0.0470 

1.53E-04 

2 

0.0030 

0.0020 

2. 066F.+05 

0.0478 

1.88E-04 

2 

0.0030 

0.0030 

2. 112E+05 

0.0490 

1.80E-04 

2 

0.0020 

0.0050 

2 . 147E+05 

0.0493 

1.58E-04 

2 

0.0030 

0.0050 

2 . 189E+05 

0.0506 

1.83E-04 

2 

0.0040 

0.0050 

2 . 234E+05 

0.0519 

2.08E-04 

2 

0.0050 

0.0050 

2 . 276E+05 

0.0519 

1.88E-04 

2 

0.0050 

0.0060 

2 . 320E+05 

0.0530 

1.64E-04 

2 

0.0060 

0.0060 

2. 368E+05 

0.0544 

1.81E-04 

2 

0.0070 

0.0060 

2 . 402E+05 

0.0557 

1.64E-04 

2 

0.0070 

0.0070 

2. 448E+05 

0.0562 

1.61E-04 

2 

0.0080 

0.0070 

2. 496E+05 

0.0576 

1.61E-04 
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TABLE  B-7.  THE  EXPERIMENT  OF  FRASIER  (1968) 


a  =  3. 

153  cm 

C  =  9. 

030  cm 

Central  Rod:  d  /a 

2  =  .023 

V  = 

2 

0.01  cm  /s  , 

p  =  0. 

,818  g/cm3 

Main 

Ring 

Increment  Rings 

Top  Bottom 

L  (g  cm2) 

tr 

"TI 

2 

0. 0000 

0 . 0000 

1 . 830E+05 

0. 0507 

2. 349E-04 

2 

0 . 0  0  1  0 

0. 0000 

1 . 878E  +  05 

0.0514 

2. 998E-04 

2 

0. 0010 

0 . 0020 

1 . 973E+05 

0. 0528 

3. 533E-04 

2 

0 . 0020 

0. 0020 

2. 018E+05 

0.  0536 

4 . 60  3  E  ~  0 4 

2 

0 . 0030 

0 . 0020 

2. O66E+05 

0.0545 

5. 138E-04 

4— 

0 . 0  0  2  5 

0 . 0035 

2. 099E+05 

0.  0552 

5. 615E-04 

2 

0 . 0  0  3  5 

0. 0025 

2. 099E+05 

0.0553 

6. 398E-04 

o 

0 . 0035 

0 . 0035 

2. 147E+05 

0.0560 

5. 672E-04 

2 

0 . 0035 

0. 0035 

2. 147E+05 

0.0562 

5 . 94  0  E - 0  4 

o 

0. 0035 

6.0045 

2. 193E+05 

0. 0570 

4 . 966E-04 

71 

0 . 0035 

0. 0045 

2. 193E+05 

0. 0570 

6 . 3  9  8  E  -  0  4 

2 

0 . 0  0  4  5 

0. 0045 

2. 239E+05 

0. 0576 

4. 584E-04 

2 

0 . 0045 

0. 0045 

2. 239E+05 

0 . 05 1*’  8 

6 . 0  92E-04 

2 

0 .005 0 

0 . 0050 

2. 276E+05 

0. 0580 

3. 285E-04 

o 

0 . 0060 

0. 0060 

2. 368E+05 

0.86  0 6 

2. 4 6 4 E - 0 4 

ii! 

0 . 0070 

0. 0070 

2. 448E+05 

0. 0623 

1 .  1  4  6  E  -  0  4 

2 

0. 0075 

0 . 0085 

2. 537E+05 

0. 0640 

8. 212E-05 

v  =3 

0.03  cm  "Vs 

>  p  = 

0.900  g/cm3 

Main 

Ring 

Increment  Rings 

Top  Bottom 

L(g  cm2) 

tr 

"TI 

0 . 0000 

0. 0000 

1 . S30E+05 

0.05O5 

2. 31  IE-04 

2 

0 . 0000 

0. 0010 

1 .  S.-8E  +  05 

0.0513 

3. 075E-04 

2. 

0 . 0  0  2  0 

0 . 0020 

2. 01SE+05 

0. 0540 

3. 915E-04 

2 

0 . 0  0  2  5 

0. 0025 

2. 051E+05 

0.0547 

4. 183E-04 

2 

0 . 0035 

0. 0025 

2. 099E+O5 

0.0557 

3. 591E-04 

2 

0 . 0035 

0. 0035 

2. 147E+05 

0.0563 

4 . 335E-04 

>*r- 

0 . 0035 

0. 0045 

2. 193E+05 

0. 0570 

4. 163E-04 

2 

0 . 0045 

0 . 0045 

2. 239E+05 

0.0578 

4.221E-04 

2 

0 . 0045 

0 . 0055 

2. 287E+05 

0. 0587 

4. 5S4E-04 

2  '  ■ 

0 . 0050 

0 . 0060 

2. 32OE+05 

0.0596 

3.320E-04 

2 

0 . 0060 

0 . 0060 

•  2. 368E  +  05 

0. 0603 

3 . 7  62E-04 

2 

0. 0060 

0. 0070 

2. 402E+05 

0. 0613 

3. O56E-04 

s-, 

0 . 0070 

0 . 0070 

2. 448E+05 

0.0622 

2 . 502E-04 

2 

0 . 0070  '• 

0 . 0075 

2. 473E+05 

0. 0630 

1.872E-04 
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TABLE  B-7. 


Increment 

Top 

0 . 0000 
0 . 0  0  0  0 
0 . 0  0  2  0 
0 . 0035 
0 . 0035 
0 . 0050 
0 . 0  0  5  0 
0 . 0  0  7  0 
0 . 0075 
0 . 0035 
0 . 0035 
0 . 0  1  05 


THE  EXPERIMENT  OF  FRASIER  (1963)  (continued) 
0.13  cm^/s  ,  p  =  0.940  g/cm^ 


Rings  2 

Bottom  L(g  cm  )  TR  TI 

0.0000  1 . 830E+0 


0. 0020 

1 . 925E  +  05 

0 . 0020 

2. 01SE+05 

0 . 0  0  2  5 

2. O99E+05 

0. 0045 

2. 193E+05 

0 . 0  0  5  0 

2 . 2  7  6  E  +  0  5 

0.0070  • 

2. 365E+05 

0. 0070 

2. 44SE+05 

0 . 0  0  3  5 

2. 537E+05 

0 . 0035 

2. 535E+05 

0 . 0  0  9  5 

41 .  6  6  3  E  +  0  5 

0. 0095 

2.71 1E+05 

0. 0504 

1 . 

280E-04 

0.0521 

1 . 

261E-04 

0 . 054  0 

1  . 

910E-04 

0.0564 

£  , 

0  2  4  E  —  0  4 

0.05  r  3 

•~i 
k_  • 

502E-O4 

0  a  0  5  O  3 

•“) 

540E-04 

0 . 06  0  9 

o 

k_  ■ 

655E-04 

0 . 0624 

£  , 

292E-04 

0. 8640 

£  , 

254E-04 

0. 0643 

349E-04 

0 . 066  3 

i ! 

833E-04 

0 . 0671 

i . 

833E-04 
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USER  EVALUATION  OF  REPORT 


Please  take  a  few  minutes  to  answer  the  questions  below;  tear  out 
this  sheet  and  return  it  to  Director,  US  Army  Ballistic  Research 
Laboratory,  ARRADCOM,  ATTN:  DRDAR-TSB,  Aberdeen  Proving  Ground, 
Maryland  21005.  Your  comments  will  provide  us  with  information 
for  improving  future  reports. 

1 .  BRL  Report  Number 

2.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related 
project,  or  other  area  of  interest  for  which  report  will  be  used.) 


3.  How,  specifically,  is  the  report  being  used?  (Information 
source,  design  data  or  procedure,  management  procedure,  source  of 
ideas,  etc.) _ 


4.  Has  the  information  in  this  report  led  to  any  quantitative 
savings  as  far  as  man-hours/ contract  dollars  saved,  operating  costs 
avoided,  efficiencies  achieved,  etc.?  If  so,  please  elaborate. 


5.  General  Comments  (Indicate  what  you  think  should  be  changed  to 
make  this  report  and  future  reports  of  this  type  more  responsive 
to  your  needs,  more  usable,  improve  readability,  etc.) 


6.  If  you  would  like  to  be  contacted  by  the  personnel  who  prepared 
this  report  to  raise  specific  questions  or  discuss  the  topic, 
please  fill  in  the  following  information. 

Name : 


Telephone  Number: 
Organization  Address: 


